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A convenient and general decomposition of the photon self-energy in a magnetized, but otherwise 
isotropic, medium is given in terms of the minimal set of tensors consistent with the transversality 
condition. As we show, the self-energy in such a medium is completely parametrized in terms 
of nine independent form factors, and they reduce to three in the long wavelength limit. We 
' consider in detail an electron gas with a background magnetic field, and using finite temperature 

field theory methods, we obtain the one-loop formulas for the form factors, which are exact to all 
orders in the magnetic field. Explicit results are derived for a variety of physical conditions. In 
the appropriate limits, we recover the well-known semi-classical results for the photon dispersion 
■^T^ relations and the Faraday effect. In more general cases, where the semi-classical treatment or the 

linear approximation (weak field limit) are not applicable, our formulas provide a consistent and 
systematic way for computing the self-energy form factors and, from them, the photon dispersion 
relations. 
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I. INTRODUCTION 

The propagation of photons in a material medium in the presence of a magnetic field gives rise to many observable 
effects. For example, in astrophysical contexts, the fact that the photons with different polarizations have different 



dispersion relations (birefringence), leads to the Faraday rotation effects that have been observed for various astro- 
physical objects (l], ||, ||, |[ || |(|. Recently 0], the observation of birefringence effects have also been reported in an 
""■-^ ' experiment that studied the millisecond pulsations of solar radio emission. In a different context, in another recent 
I laboratory experiment]^], an effect that mimics the Faraday rotation has been found. In this case, using two input 
^p- beams (a control and a probe) it was shown that a Faraday-like effect is induced, for both circular and linear po- 
O i| larizations of the control beam. From a theoretical point of view, the common theme in all these situations is the 
propagation of a photon in matter in the presence of an external field, being a constant magnetic field in one case, or 
a plane wave electromagnetic field in the other. 

The subject of the propagation of photons in magnetized plasmas is of course a well studied one, and it is amply 
covered in standard astrophysics and plasma physics textbooks]^, |l^|. The traditional approach to this subject has 
been based on the equations of classical kinetic theory and related semiclassical methods. On the other hand, it is 
■ commonly accepted that the field-theoretical methods are generally useful in this kind of problem, since they are 
applicable to a wider range of physical situations for which the semiclassical methods breakdown. Such situations 
include, for example, matter under extreme conditions, such as highly relativistic and/or highly degenerate. 

For the weak field limit a field-theoretical treatment has been given in Ref . jll| . The assumption there is that the 
photon self-energy can be expanded as a power series in the external magnetic field and that the linear terms are 
the relevant one. There exist physical situations that lie outside of the scope of such approximation and require a 
more general treatment. For example, the arguments for the existence of super critical fields (> 4.4 x 10 13 G) in 
nature have grown stronger recently with the observations of few Soft Gamma Rays Repeaters (SGR) and Anomalous 
X-Rays Pulsars (AXP), which are very likely magnetars[p^| [l3| , that is, isolated neutron stars with surface magnetic 
fields of order 10 14 — 10 15 G. A model for extragalactic gamma rays burst in terms of merger of massive binary stars 
suggests also that magnetic fields up to the order of 10 17 G|l4| may exist. In the context of the Early Universe, 
very large magnetic fields (10 23 G) may be generated |l5) during the electroweak phase transition due to gradients in 
the Higgs field. Thus, there are environments of interest - that involve, in addition to matter in extreme relativistic 
and/or degenerate conditions, strong magnetic fields - for which neither the semiclassical methods nor the weak field 
approach are directly applicable. 

Our purpose here is to have another look at this subject, keeping the above points in mind. To this end, we 
consider the photon propagation in an electron gas in the presence of a constant external magnetic field. We derive 
the general decomposition of the photon sclf-cncrgy in such a medium, in terms of the minimal set of tensors consistent 
with isotropy and the transversality condition. We show that the self-energy is completely parametrized in terms of 
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nine independent form factors, and that in the long wavelength limit only three are independent, with the rest 
being expressed in terms of them. Using the (real-time) finite temperature field theory method, we obtain the one- 
loop expressions for the form factors, expressed as integrals over the distribution functions of the particles in the 
background. Those formulas are valid to all orders in the magnetic field and are explicitly evaluated for a variety 
of conditions of physical interest. As an application we determine the photon dispersion relations and calculate the 
Faraday rotation for plane polarized light, for various cases. They reproduce the well-known semi-classical results 
when the appropriate limits are taken, but remain valid for more general situations. 

In Section II we give the general decomposition of the photon self-energy in terms of the nine independent form 
factors. We also collect there several kinematical relations that are useful in later stages of the calculations. The 
one-loop formulas for the photon polarization tensor are derived in Section III. These are used in Section [fy| to 



discuss the dispersion relations, focusing on the long wavelength limit as a special case. The explicit formulas for the 
independent form factors are given in Section ^ for different possible conditions of the background electron gas, and 
for various regimes of interest including the low frequency regime and the weak-field (linear) limit. Our conclusions are 
summarized in Section VI , while in Appendix ^ we summarize the conventions that we use regarding the Schwinger 



formula for the electron propagator in an external magnetic field, and Appendix |B| contains the derivation of an 
integral formula used in the calculation. Some of the details of the derivation of the low-frequency formulas, and the 
weak-field formulas, for the self-energy form factors are shown in appendices ^ and respectively. 

II. KINEMATICS 

A. General decomposition of the polarization tensor 

In the absence of the magnetic field, the photon self-energy tt^ depends in general on the photon momentum q^, 
and on the velocity four- vector of the medium u^. In the frame of reference in which the medium is at rest, w M has 
components given by 

u"=(l,6), (2.1) 

and in that frame we write 

g" = (w,<3). (2.2) 

In the presence of a magnetic field, but otherwise an isotropic medium, ir^ depends in addition on the vector 6 M that 
is determined by the magnetic field. The vector b^ is defined such that, in the frame in which the medium is at rest, 

6" = (0,6), (2.3) 

where we denote the magnetic field vector by 

B = Bb . (2.4) 

For a given photon momentum vector we define the unit vectors (i = 1, 2, 3) by writing 

Q = Qe 3 , (2.5) 

with ei t 2 chosen such that 

ei • e 3 = ei • e 3 = , 

e 2 = e 3 x ei . (2.6) 

In addition, for the problem that we are considering in the present work, without loss of generality, we can choose the 
vectors ei,2 such that b lies in the 1, 3 plane. Thus we can write, 

b = cos #£3 + sinflei , (2.7) 

where 



cos 9 = Q ■ b . 



(2.8) 
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We now introduce the vectors 

X. 



n _ n (u ■ q)q t ' 



*7 j 



X% = e^q a bpu~, 
X? = e^X 2a q 0Ul , (2.9) 



which satisfy 



q-Xi = 0, 

XfXj = (i^j). (2.10) 



Therefore, the X? form a basis of vectors orthogonal to g M , and the nine bilinear combinations XjfXj form a complete 
set of tensors in terms of which the photon self-energy can be decomposed. To exploit this more fully, it is useful to 
define the (normalized) vectors 



which satisfy 



«f ( « = 7=rx' (2 - n) 



u(Q)-q = o, 

= - E ^i(Q>uj(Q) , (2.12) 



q- 



i=l,3 



and which take the form 



e$*(Q) = (0,e 



4(Q) = (0,e 2 ), 

e£(Q) = -L=(Q, w g 3 ) (2.13) 



in the rest frame of the medium. 

Whence, in the most general case, the photon self-energy can be expressed as a linear combination of the bilinears 
tfe'i in the form 



involving the nine independent coefficient functions n^' (u>,Q). Furthermore, as will be seen in Section IVB, in the 



long wavelength limit all the (ui , Q) can be determined in terms of just three independent functions [i.e., Eq. 



(4.33)]. This fact will be the starting point in Section IV to obtain the polarization vectors and dispersion relations 



of the normal modes. 

B. Additional kinematic relations 

It will prove to be useful to introduce the following tensors that are transverse to u^, 

Qui; bfj,bij , 

P„ v = i^v a pb a uP . (2.15) 

They have the orthogonality properties 

QR = QP = 0, (2.16) 
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and they satisfy the multiplication rules 



-,2 



R? — R , 
P 2 = R, 



as well as the normalization conditions 



In addition, we define 



RP = P, (2.17) 



Q% = 1, 

i?% = 2. (2.18) 



= u^K - b^Uy , (2.19) 



which satisfies 



UR = RU, 
U» x U Xv = u»u v - b"b v , 

= <f„ - R\ , 

\p^€, luaP = iU a/3 . (2.20) 

In any case, any four-vector u M can be decomposed according to 

w" = ^ + <, (2.21) 

where 

< = R» u v" = (0,i? ± ) , 

w|| - w"-< = (u°,U||), (2.22) 

with and iT|| being the three-dimensional components of v that are perpendicular and parallel to b, respectively. 
In particular, for any two such vectors, 

v± ■ v'j_ — v ■ v' — (v ■ u){v' ■ u) + (v ■ b){v' ■ b) , 

v v v[ = (v ■ u)(v' ■ u) - (v ■ b)(v' ■ b) . (2.23) 

III. PHOTON SELF-ENERGY 

In one-loop, the 11 element of the photon thermal self-energy is given by 

j^Tt [ lti iS e (p + qhuiS e (p)} , (3.1) 

where S e stands for the 11 component of the electron thermal propagator. It can be decomposed in the form 

S e (p) = S Fe (p)+S Te (p), (3.2) 

where Sp e is the propagator in vacuum in thepresence of the magnetic field while St c incorporates the effects of the 
thermal background. As shown in Appendix ^, Sp(p) can be written in the form 



iSp{p) = / drG(p,s)e 
Jo 



ir<&(p,s) — e 



(3.3) 
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where 



with 



and 



*(p, s) = p 2 -m 2 e + (s- 1 tan(s) - l) R^jf , 
G(p,s) = \f 1 (p,s)+i-f 5 y 2 (p,s)+m, e + m e tan(s)ij 5 '4fi 

Vi M (p, s) = sec 2 (s)p M + tan 2 (s) [(p • b)b^ - (p ■ u)tt M ] , 
V2 M (P> S ) = tan(s) [(p • fyup - (p ■ u)6 M ] 

s = lelBr. 



Using the relation sec (s) = 1 + tan (s), Eq. ( J3 . 5| ) can be rewritten as 

Vip(p,s) = [g^u + tan 2 {s)R fll/ ] p v , 
V2|*(p, s) = tan(s)[/ M , y p' y . 



o<T*(Pi»)- 



On the other hand, the thermal part is given by 

/oo 
drG(p, s)e l 
-oo 

where 

Ve(p) = 6(P ■ u)f e (p ■ u) + 6{-p ■ U)fe(-P ■ u) 

with 

fe(x) 
fe(x) 



(3.4) 

(3.5) 
(3.6) 

(3.7) 
(3.8) 
(3.9) 



1 

e /3(x-/i e ) _)_ l 
1 

g/3(.T + ^ e ) _|_ -y 



(3.10) 



Here /? stand s for the inverse temperature and /i e the electron chemical potential. 

When Eq. (|3.2| ) is substituted into Eq. (|3.l|), we obtain four different terms. Only those two that contain one factor 

of SVe and Sxe each, contribute to the real part of the self-energy, which we denote by Tr^J^\ and is the quantity in 
which we are interested. Thus 



d 4 p 

(2tt)* 



(p) / dre x{ ^ / dT'e xip+,1 ^TT[^G(p + q,T')^G(P,r)] 



+ Ve(p + q) dre x ^ rfr'e^+^Tr [ 7 „G(p + g, r')7^(p, • 



where, to simplify the notation, we have defined 

A(p,r) = ir$(p, s) - e|r| 



(3.11) 



(3.12) 



In the term involving the factor r] e (p + g), we make the change of variables p — > p — q, and r <-» r'. Then, using the 
cyclic property of the trace, we obtain 

4J f \u>,$) = Ue 2 J ^Ve(p)J dre x ^ |jf dr' [e^+^U^p, g) + e^-^U^p, -g)] | , (3.13) 



where 



g) = ^ Tr [7m G (p + ?! T 'hvG{p, t)] 



Using the expression for G given in Eq. (3.4), by straightforward evaluation of the trace we obtain 

L — TS S ^ + 



(3.14) 



(3.15) 
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where 



LW (p, 9 ) = + V{pVi v - Vi ■ v{ 9tw - v^v{ v - v 2 ^v 2v + u 2 • u 2 V , 

+ m 2 ^,, + to 2 tan(s) tan(s') (g^ - 2u il u v + 2b fl b 1/ ) , 
(P, q) = enu a pV{ a v£ + e^pV^Vf + m 2 [tan(s) - tan(s')] e^ a pu a b^ . 



(3.16) 



with the convention that e 0123 = i ; and using the notation s' = |e|Br', — Vi fl (p + q,s'), and similarly for V 2fi . 
Using Eq. (3/7), L^li can be expressed in the form 

(P, l) = 2m l tan(s) tan(s')^V + m l I 1 _ tan(s) tan(s')] g M „ - g^ [1 + tan(s) tan(s')] (p 2 + p ■ q) 
- g^ u [tan 2 (s) + tan 2 (s') - tan(s) tan(s') + tan 2 (s) tan 2 (s')] (p^ + pj_ ■ q±) 
+ { [g^a + tan 2 (s)i? MQ ] [g v p + tan 2 (s')^/3] - tan(s) tan(s')^^/3 + (A* ^ f)}p a (P + ?) /3 , 
lW(p, q) = e^p [ta,n(s)g aX U 0p + tan(s')£/ Q V P + tan(s) tan 2 (s')i? QA C/ /3p + tan(s') tan 2 (s)U aX R 0p ] Px (p + q) p 
+ m 2 (tan(s) - tan(s')) W^"^ , (3.17) 
where p^ and g p are defined according to Eq. (2.22). It is seen from Eq. (3.17) that 

4r ) (-P,«) = 4 s / ) (p,-3)- (3.i8) 



Using this and the fact that A(— p, s) = A(p, s), Eq. ( 3.13 ) can be written in the form 

d 4 p 



/p-u>0 (271") 

s^lWfp, g)[/ e (p • u) - hip ■ «)] - (g -> -g)] } 
The integrals over r, r' can be expressed in the form 



dre x( ^ / dr' { [e A ^ T UW(p, g)[/ e (p • u) + f s (p ■ «)] + (g -> -g) 



(3.19) 



dre^^ / dr'e^+^L^ip, q) = -2m 2 #i(p)Ji(p + g)i^„ + m 2 ^ [X (p) Jo(p + ?) + #i(p)Ji(p + g)] 
Jo 

- g MV [#„ (p) Jo (p + g) - K x (p) Ji (p + g)] (p 2 + p • g) 

[-if 2 (p) J (p + g) - K {p)J 2 {p + q) 
+K 1 {p)J 1 (p + q) + K 2 (p)J 2 (p + q)] [p\ + p x ■ g x ) 

+ [gna9upKo(p)J (p + g) - R^ a gvpK 2 (p)Jo(p + g) 
-g^ a R v pK {p)J 2 (p + q) + R fia Ru/3K2(p)J2{p + g) 

+ U^U^K^ppiip + q) + (/! «-> i/)] p<*(p + g)' 3 , 

[^i(p)Jo(p + g)g aA ^ p + k (p)Mp + g)c/ Q V p 

- Xi(p)J 2 (p + q)R aX U f3p - K 2 {p)Ji(p + q)U aX R^\ p A (p + g) p 
+ m'te^u ^ [Xi(p) J (p + g) - K {p).h{p + q)} , (3.20) 



dre A(p ' T) / dr'e A(p+9 ' T ' ) 4 a j(p,(7) - 



where 



with 



Jn(p) 



-»)* 



dse^e- io " taI1 Wtan n (s), 



(3.21) 



9 -^9 

= -pi _ p[_ 
P ~ \e\B \e\B ' 



\B ( 



|e|S 



P\\ 



m„ +ie) = 



\e\B 



2 ^2 2 i 

p -py -m e +ze 



(3.22) 



The K n are defined by the same integrals as in Eq. ( 3.21 ), but with the limits of integration being — oo < s < oo 
instead. Therefore, the following relations are easily obtained 



K n (p) = J n (p) + J*(p) 



(3.23) 



Furthermore, from Eq. (3.21) it is easily seen that 



Jn+l(p) = 



dJnjp) 

da r . 



(3.24) 



and therefore, only Jo needs to be evaluated explicitly. As shown in Appendix ||, it can be expressed in the form 

D e (a p ) 



MP) 



E 



\e\B ^ V^p - 2£ 

where the D n (a p ) are given in terms of the Laguerre polynomials by 

D n (a) = {-l) n er a [L n (2a) - L„_i(2o)] 

with the convention 

L_i = 0. 



From the relations in Eq. (3.24) we then obtain 



Hp) 



D'Mp) 



i D"(a p ) 

\e\B ^ i> v -2e' 



where 



dD n (a) 

da 
d 2 D n (a) 

da 2 



Furthermore, from Eq. (3.23), it follows that 

oo 

K (p) = 2TrY,D n (a p )S(pj~ml-2n\e\B). 

71=0 
OO 

K^p) = 2TTY / D' n (a P )S(pj~m 2 e -2n\e\B). 

oo 

K 2 (p) = 27rJ2KM^P\\-m 2 e -2n\e\B). 



Substituting Eq. (3.20) into (3.19), and using Eq. (3.30) wc then obtain 

d 3 p 



{ [MW(p, q)[f e {E n ) + f s (E n )} + (q^ -q) 



"J {2ir) 3 2E n 
+ [N$(p, q)[f e (E n ) - f s (E n )\ -{q^ -q)] } 



where 



with 



P° = E n 



(3.25) 

(3.26) 
(3.27) 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



E n = J p.? +ml + 2n\e\B , 



(3.33) 
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and 



M$(p, q) = -2m 2 e D' n (ap)Ji(p + q)R„ v + m\ [D n (a p )J (p + q) + D' n {a p )J x {jp + q)] 

- [D n {a p )J (p + q) - D' n (a p )Ji(p + q)} (p 2 + p ■ q) 

- 9»» [- D n( a p)Jo(p + q)~ D n (a p )J 2 (p + q) 
+D' n {a p ) Jt{p + q) + D' r [(a p )J 2 (p + q)] (p 2 ± + p ± ■ q ± ) 
+ [9na9 v pD n (a p )J {p + q)- R m 9u^D'^{a p )J Q (j) + q) 
-9^ccRupD n {a p )J2(p + q) + R^aRvpD'^dp) J 2 {p + q) 
+ U t , a U 1J pD' n {a p )J 1 (p + q) + (fi «-► v)]p a {p + qf , 

N$(p,q) = ie^[D , n (a p )J (p + q)g aX U^ + D n (a p )J 1 (p + q)U aX g^ 

- D' n {a p )J 2 (p + q)R aX U^ - D^Mp + q)U aX R^) p x (p + q) p 
+ mlie^ a f 3 u a b f3 [D' n (a p )J (p + q) - D n (a p )J\{p + q)] . 



(3.34) 



Using Eq. ( 3.31 ) as the starting point, the form factors that enter in the dispersion relations of the propagating modes 
can be computed, as we show next. 



IV. DISPERSION RELATIONS OF THE NORMAL MODES 



A. General case 



As already argued in Section II, the photon self-energy can be expressed as shown in Eq. (2.14). The nine indepen- 
dent coefficients ^'^(u), Q) are then determined to one-loop order by 



(4.1) 



With TT$ f) (uJ,Q) given by Eq. (|3.31|). While e^(Q) and £%(Q) are independent of Q, e^(Q) is not. We will denote by 



(0) the basis vectors for Q = 0, and therefore 



< a (0) = <a(Q) 
e£(0) = (0,e 3 ). 



It is useful to note that 



and by using 



we can also write 



e£(Q) = 4f ^ 

v q 



4(Q) 



<t = wu" + Qe£(0) , 
Q u 



4(0) 



The polarization vectors of the various propagating modes can be expanded in the form 

e= £ Oi(3K(3), 

i=l,3 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



where the coefficients ai(Q), and the corresponding dispersion relations, are determined by solving the equation 



In matrix notation, Eq. (4.7) can be written in the form 

(q 2 -n)a = 0, 



(4.8) 
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where 



a 




(4.9) 



and II is the matrix formed by the coefficients n^K The formulas for the functions are obtained by applying 
Eq. (4.1). For that purpose, it is useful to note that we can make the replacement 



4(Q) 



q e£(0), 



(4.10) 



which follows by using Eq. (4.5) and the transversality condition satisfied by irffi^ . Thus we obtain 





_ % {eff) o&( W) Q) 


{a3) (u,Q) = 




i3a) (u,Q) = 




LU 


{33) (co,Q) = 


_4 7r (e//)33 (Wj Q ); 



(4.11) 



wher e the indices a, b can take the values between 1 and 2. While Eq. (4.11) allows us to read the elements tt^' off 
Eq. (3.19) by inspection, the problem of finding the dispersion relations in the general case is a formidable one. For 
this reason we now consider in some detail the so called long wavelength limit, which is a particularly important case 
that holds in a variety of physical applications. 



The point to stress here is that the nine elements n^' determined by Eq. ( 4.11 ), form a complete set of independent 
functions that parametrize the photon self-energy in the most general way that is consistent with the transversality 
condition. In particular, this characteristic remains valid independently of any appro xima tion that may be used to 
compute the elements ir^fi^i that must be substituted in the right-hand side in Eq. (4.11). 



B. Long wavelength limit 

The long wavelength (Q — > 0) limit is particularly important and we consider it separately as a special case. 
We stress that, while we are considering this particular limiting case, we do not make any assumptions about the 
conditions of the gas or the magnitude of the magnetic field. Therefore, the formulas that we obtain below hold for 
relativistic or non-relativistic gases, whether they are degenerate or not, and for any value of the magnetic field. This 
limit is valid under the condition 

uj^>v e Q, (4.12) 

where v e stands for the average velocity of an electron in the gas. In the limit Q — > 0, we can write 

qf* = lou* 1 , (4.13) 

and therefore only and are independent vectors. The most general form of it^J consistent with the transver- 
sality condition in this case is then 

7r£,'fl (w, Q -» 0) = 7r T (c;, Q -» 0)R^ + tt l (lu, Q -> 0)Q^ + tt p (lu, Q - 0)P M „ , (4.14) 

where R^v, P^v and Q^ v are defined in Eq. ( 2.15| ), and the functions ttt,l,p are determined from the one- loop expression 
for TTffJ by means of the projection formulas 

n P (Lu,Q) = -l P ^(fiff){ u ,Q). (4.15) 
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Although we are not in dicati ng it explicitly, the functions ttt,l.p(u,Q) depend, in addition, on the magnetic field B. 
Substituting Eq. ( 3.19 ) into ( 4.15 ) we obtain 



tt t ,l(w,Q) = 4ie 2 ^ 



<J " 1 ' [f e (E n )+f- e (E n )} \4 n l(p, q ) + ( q ^-q) 



{2n) 3 2E n 



where 



{f e (E n ) - f e (E n )] \lP (p, g) - (g - -q) 



{2n) 3 2E [ 



4 n) ee \r^m$ , 



(4.16) 



r(») _ 



r u M (n) j 



r (n) 



,_pnvja(n) 

2 V ' y 



By direct calculation, it is straightforward to obtain, for any g M , the following expressions 
4 n) = [£> n (a p )J (p + g) - D' n (a p )Ji(p + q)]h T (p, q) , 



(4.17) 



where 



r(") 



[D n (a p )J (p + q) + D"(a p )J 2 (p + q) - D n (a p )J 2 (p + q) - £>"(a P ) J Q (p + q)]h L1 (p, q) 



- [D n (a p )J (p + q)+ D' n (a p )Ji(p + q)]h L2 (p, q) , 
j£° = [£>„(a p )Ji(p + g)-£);(a p )Jo(p + g)]ftTfeg) 1 

/it = (p ■ b) 2 ~ (p ■ u) 2 + (p ■ b)(q ■ b) - (p ■ u)(q ■ u) + m 
hm = P± +P± ■ gj_ , 



(4.18) 



2 

e ' 



>L2 



(p • u) 2 + (p • 6) 2 + (p • 6)(g • b) + (p • u)(g • u) - m\ 



(4.19) 



In order to evaluate the functions ttt,l,p in the Q — > limit, as indicated in Eq. (4.14), we proceed as follows. 
Remembering that p° = E n , in this limit 

ii p+q = -^{lu 2 + 2uE n + 2n\e\B) 7 (4.20) 



and then from Eq. fl3.25| ) 



14 s 



Jo{p + q) = 



Dn(a p ) 



2ujE n + 2{n- l)\e\B) 



(4.21) 



with similar formulas for J\ t2 {jp + g), but with Di(a p ) being replaced by _D^(a p ) and Z?"(o p ), respectively. Taking the 



Q — * limit of the expressions in Eq. ( 4. IS ) and substituting them in Eqs. ( 4. IS ) we then obtain 



r (n) 



rW 



(D n {a p )D e (a p ) - D' n {a p ) D' e (a p )) (coE n + 2n\e\B) 



ia p \e\B 



lo 2 + 2uE n + 2(n- £)\e\B 

D n {a p )D e (a p ) + £%(a p )£>?(a p ) - £> n (n p )£>?(n p ) - D^(a p )D e (a p ) 
uj 2 + 2ujE n + 2{n- l)\e\B 



An) 

1 p 



[D n (a p )D e (a p ) + D' n (a p )D' e (a p )} [2E 2 n + uE n - 2(m\ + n\e\B)} 
lu 2 + 2uE n + 2(n - £)\e\B 

{D n {a p )D' £ (a p ) - D' n (a p )D e (a p )) (uE n + 2n\e\B) 
u> 2 + 2uE„ +2{n- £)\e\B 



(4.22) 
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With the help of the identities, 



^ (L n — L n _i) — —L n _i , 
ax 



the following is easily derived 



a(D n (a)-D'^(a)) =2nD n (a), 



(4.23) 
(4.24) 



so that the formula for / 



(«)l 



r(n) 



becomes 



dp 



4n£D n (a p )D e (a p ) 



uj 2 + 2ujE n + 2(n - £)\e\B 
D n (a p )D e (a p ) + D' n {a p )D' e (a p )} [2E 2 n + wE n - 2(m\ + n|e|B)] 



; 2 + 2w£;„ + 2(7i-^)|e|B 



The integrand in Eq. ( 4.16 ) depends on only through a p , so that we can replace d 3 p — > dpii (7r|e|-Bda p ) 
the help of the normalization condition satisfied by the Laguerre polynomials 

dxe x L m (x)L n (x) = S n>m , 

where 5 nm is the Kronecker symbol with the convention 

3n,n = for n < , 

the following hold 



daD n (a)D m (a) 



daD' n (a)D' m (a) 



(-1) 



n-\-rn 



[^n,m ^n— l.m $n,m— 1 H~ ^n— 1 



(-1) 



n+m 



[^n,m ~h ^n— l.rn ^n,m— 1 ~t~ <^n— l,m— l] 



daD n (a)D f m (a) = - 



(-1) 



[^n,m ^n — l,m ~t~^n,m — 1 ^n— l,m— l] 



da 



D n (a)D m (a) 



Thus, defining 



70c (w) = 7Tx (w, <2 — > 0) (for X = T, L, P) , 



by substituting Eq. (4.22) in Eq. (4.16) and using the above integration formulas we obtain 



7T X (w) 



\ e ? B \r Mr a 



where 

(0)/ 
7Ty (W 



(n/Q) / 



(0)/ 
7Tp (w 



(n#0) , 
7Tp (W 



(0)/ 



(n^O) , 



dp||(/e(^)+/e(£ ))- 



w 2 + 2w£ - 2\e\B 
f 'p-ifeiEn) + fs{E n )) {uE n + 2n\e\B) 

dp\\(fe(Eo) — fe(Eo)) 

dp 



E n 



uj 2 + 2ujE -2\e\B 
{fe(En) - fe(E n )) (ujE n + 2n\e\B) 



n=0 

+ (uj -> -U)) , 
1 

uj 2 + 2ujE n - 2\e\B + uj 2 + 2uE n + 2\e\B 
- (uj -> -w) , 

1 1 



1 



w 2 + 2u)E n - 2\e\B uj 2 + 2uE n + 2\e\B 



^(fe(Eo) + MEo)) + (a; - -w) 



r ^i(f e ( En )+ME n )) 



8n\e\B + 4ml 
uj 2 + 2uE n 



+ (uj — > — w) 



(4.25) 
With 

(4.26) 
(4.27) 



(4.28) 
(4.29) 
(4.30) 



+ (uj — > -cj) 



(w — > — oj) 



(4.31) 
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Since we are concerned only with the real part of the self-energy in this work, the above integrals are actually defined 
in the sense of their principal value part. 

As we discuss next, the formulas in Eq. ( [4.31 ) allow us to study the dispersion relations in the long wavelength 
limit. In this case Eq. (4.8) becomes 



[q 2 - n ]a = . 



(4.32) 



Here the matrix n is formed by the coefficients 7rW' that are determined from Eq. (4.11), but using the formula for 



r (e//), 



(uj, Q — > 0) that is given by Eq. (4.14). Thus, using the kinematical relations given in Sec. II, by straightforward 



algebra we find 



n n = 



/ 7Tt(w) cos 2 9 + ttl(lu) sin 2 9 — iirp(uj) cos 9 [n l (lu) ~ ttt{uj)} sin 9 cos 9 ^ 



iTTp(uj)cos9 7r<r(o;) — in p{lo) s\i\ 9 

\ ^f-[-K L {uj) - tt t {uj)} sin#cos# ^f-iir P (uj) sin9 ^ [tt t (lu) sin 2 9 + tt l (uj) cos 2 9] ) 



(4.33) 



It should be remembered that in Eq. ( 4.32 ), a depends on Q. In what follows, we consider some particular cases. 



1. Q = 

We consider first the case the zero momentum limit of the dispersion relations, which corresponds to a photon 
traveling with a vanishingly small momentum in an arbitrary direction. The matrix n in this case is given as in Eq. 
(4.33) but with the replacement q 2 /uj 2 — ► 1. The eigenvectors of the matrix in this (Q — 0) case are found to be given 
by 



a± 



OiL 




with corresponding eigenvalues 7r^ ± tt p and 7r^, respectively. Therefore, from Eq. (4.6), in the Q 
polarization vectors of the normal modes are given by 

££(0) = cos0e£(O)±ie£(O) -sin0e£(O), 

f£(0) = sin0e^(O)+cos0e£(O), 

and the corresponding dispersion relations are determined by solving 

UJ 2 — [TTT ± 7Tp] = , 



respectively. 



L0 2 -TT L = 0, 



2. Propagation parallel to the magnetic field 



(4.34) 
limit the 

(4.35) 
(4.36) 



This corresponds to set 







(4.37) 



in Eq. (4.33). However, in contrast to th e prev ious case, we do not set Q — (although we still assume that Qv e <C u>, 



which is the condition under which Eq. (4.33) holds). In this case 



7T T — ijTp 

n (c^Q) = ( in° P n T 

<7T° 



(4.38) 
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and therefore 



a± = 



(4.39) 



with corresponding eigenvalues 7r^ ± 7Tp and j& 7r £> respectively. The polarization vectors of the normal modes are 
given by 

e ± (Q) = i(e?(Q)±ie£(Q)) , 

- e%(Q), (4.40) 
and the corresponding dispersion relations are determined by solving 

UJ 2 - Q 2 - [7T T ±7Tp] = 0, 

uj 2 -tt l = 0, (4.41) 

respectively. One of the characteristics of this solution is the fact that the longitudinal mode dispersion relation is 
independent of Q, which is a well known feature of this case. 



3. Propagation perpendicular to the magnetic field 



In this case 



tt/2 



and, as in the previous one, we maintain Q ^= 0. Then, from Eq. (4.33), 

/7T L 



n (w,Q) 



7T T -i^ir p 



(4.42) 



(4.43) 



The dispersion relation of the longitudinal mode [with the polarization vector ££(Q) = e^{Q)] is obtained by solving 

uj 2 - Q 2 - 7T L = 0. (4.44) 

For the transverse modes it becomes 

(lu 2 -tt t ) 2 -Q 2 {lo 2 -n T )-Tr 2 P = 0, (4.45) 
which can be written in the alternate form 



Q 2 (i-5) 2 -(^) 2 

W 2 1-2$. 



(4.46) 



The corresponding eigenvectors are linear combinations of e% ^iQ)- 

For an arbitrary direction of propagation, the longitudinal and transverse polarizations mix, so that the propagating 
modes are neither purely longitudinal nor purely transverse. In this case, a perturbative solution, that is valid for 
sufficiently small Q, can obtained. 



EXPLICIT FORMULAS 



Although the formulas in Eq. ( 4.31 ) have been obtained by taking the long wavelength (Q — > 0) limit, no ap- 
proximations have been made. However, they can be calculated explicitly for a variety of situations, under certain 
approximations. Here we specifically consider the case in which u> satisfies 

w«v^B (5.1) 
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and 



w < 2(E e ) , 



(5.2) 



where (E e ) stands for an average value of the energy of an electron in the gas. Thus, if the gas is non-relativistic, the 
cond ition holds for uj <C m e . If the gas is extremely relativistic, it also holds for to > m e , subject to the condition in 
Eq. (5^). Together with the condition given in Eq. ( 4.12 ), this implies that we are considering a regime in which 



(5.3) 



to which we will refer as the low frequency regime. 



Then, under the assumption that the condition in Eq. (5.1) holds, the formulas for tt^ p given in Eq. (4.31) reduce 



to 



4 n) H = (2 - 5 nfi )iu 2 



d Pll (f e (E n ) + f s (E n ))- 



l E 2 



1 B 2 



E n 



n\e\B 

E n 



TT ( p\u) = [2-8 nfi ){u\e\B) / d n {f e {E n ) - fe{E n )) 



and similarly, 



l El 



2 B 2 



2n\e\B 
El 



7r< n) («) = (2 - S n , ) I dm (f e (E n ) + f s (E n )) 



(ml + 2n\e\B) 
IP 



(5.4) 



(5.5) 



assum ing Eq. ( p.2\) . We should note here that the judicious application of the low freque ncy regime conditions to Eq. 
( 4.31 ) requires some care. The procedure we have followed to arrive at Eqs. (5.4) and (5.5) is outlined in Appendix 



B 

We stress once more that these formulas hold for any conditions of the electron gas, whether it i s relativistic or 
non-relativistic, and degenerate or not, and for any value of the magnetic field, subject to Eq. (5.1). However, for 
specific conditions of the gas, these expressions reduce further to simpler formulas. We consider several examples 
separately. 



A. Nonrelativistic gas 



In this case the temperature and chemical potential of the electron gas are assumed to be such that 

it) 

I3m e > 1 , — > 1 . (5.6) 
fx 

This implies that 

h « , (5.7) 

while for / e , only values those values of p\\ and n for which 

E n ~ m e (5.8) 

contribute significantly in the integrals because otherwise f e becomes small. Therefore, to the lowest order in l/m e 
we can replace 



E n — > m e 



(5.9) 



in the integrands in Eqs. (5.4) and (5.5), and retaining the 0(l/m e ) terms leads to 



U0 2 — UJg ' 



u l (uj) = nl, 

TTp(Uj) = —z T , 



(5.10) 
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where 



L0 B 



\e\B 



e 2 n e 



(5.11) 
(5.12) 



with 



\e\B 
4tt 2 



OO 

E( 

n=0 



(2 -*„,„) / dp\\f e {E n ) 



(5.13) 



The formulas in Eqs. (5.11), (5.12) and (5.13) are the standard expressions for the cyclotron fre quenc y, the plasma 
frequency and total electron number density respectively. Moreover, using the results of Section IV B| , the formulas 
in Eq. ( [5.1C| ) allows us to reproduce very simply the classic results jl^ for the photon dispersion relations in the case 
under consideration. For example, if the photon propagates parallel to the magnetic field, the dispersion relations 
given by Eq. (4.41) become 



lo 2 -Q 2 



(U - LO B ) 



0. 



uj 2 -Q 2 - 



(lj + lo b ) 



= 0, 



(5.14) 



The formulas in Eq. (5.10), and whence those in Eq. (5.14), represent the leading term in powers of l/m e , and they 
neglect entirely the momentum-dependent terms. The latter can give non-negligible corrections at higher temperatures 
or chemical potential. However, while the leading terms given above depend only on the total number densities and 
not on the shape of the distribution function, the same is not true with the momentum-dependent corrections. Thus, 
for example, the corrections are different if we consider the gas to be degenerate or classical. 

In what follows we determine the corrections specifically for the case of a classical gas. Thus we can put 



P\'\ 



2m P 



11LOB, 



and for the distribution function 



fe 



-/3(£„-a0 



We consider first ttt,p- Using Eq. (5.15) we write the denominators of the integrands in Eq. (|5.4|) in the form 



lo 2 EI 



e 2 B 2 



Itiiob \ 
m e J 



(5.15) 



(5.16) 



(5.17) 



and similarly 



1 



n\e\B 
E n 

2n\e\B 
E 2 



m e 1 + 



1 



(5.18) 



up to terms that are most linear in /m 2 or ws/m e . The integrals over are carried out very simply and for the 
sums over n we use the formulas 



so = £• 



1 



n=0 



-/3nu B _ 



1 - e-P"B ' 

, — 0UJB 



n=l 



(1-e-^B) 2 4sinh 2fe 



(5.19) 



16 



Then eliminating the chemical potential in favor of n e by means of Eq. (5.13), which in the present case yields the 
relation 



e\B\ 2mn 



4tt 2 



ft 



-(250-1): 



e\B\ 2mn 



4tt 2 



ft tanh \ ' 



(5.20) 



this procedure yields 



ttp(uj) 



u 2 n 2 



(1 + S ) 



^|(1-^) 



or — UJ 



B 



— W(*o + «i), 

(uj 2 ~w 2 b ) 2 



(5.21) 



with 



5o 



2ujb 
m e sinh /3ojb 



(5.22) 



Notice that for Pub small, Si ~ 2<$o- As /Jwb becomes larger, 81 gives the O( 0u!% /m e ) correc tions to ttt (<*>). 

Turning the attention now to ttl, we expand E n in the integrand of Eq. ( |5.5| ) using Eq. ( 5.15| ). Retaining terms 
that are most linear in p,f /m 2 or us/m,, and proceeding as above with the integrals and sums over n, 



In the limit Plub — > 0, this formula reduces to 



1 - 



UJ B 



2m e ft m e sinh ftujs 



1 



2m e ft 



(5.23) 



(5.24) 



which is the standard tempera ture correction to the plasma frequency (in the non-relativistic, classical regime). 
Otherwise, the last term in Eq. ( 5.23 ) gives the corrections of 0(P 2 u>%/m e ) due to the presence of the magnetic field. 



B. Relativistic gas 



This limit corresponds to the conditions 



1/P,(jt, ^/\e\B > m e , 



so that we can effectively set 

and we specifically consider the situation in which 



Under these conditions 



1 5E n \e\B 



«1. 



E n 5n E 2 

(because E n is of order 1//3 or ji), and we can then carry out the sums over n by making the replacement 



£-> dn 



2n\e\B ' 



(5.25) 
(5.26) 
(5.27) 
(5.28) 

(5.29) 
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where we have defined 



pi = 2n\e\B. 



(5.30) 



With this substitution, the formula in Eq. (5.13) for the electron number density, and the analogous one for the 
positron, reduce to the standard expressions 



d 3 p 



fe. S (E) , 



where, using Eq. ( [5.3C| ), 



E 



Pi = \p\ 



(5.31) 



(5.32) 



factors of (3yJ\e\B and/or yJ\e\Bj\i. Similarly, substituting Eqs. fl5.4| ) and (5.5) into (4.30), 



In Eq. (5.31) we have neglected the contribution from the n = term, which is smaller than the dominant terms by 
d/or^ 

ttt (ui) 

7TL M 



2eW 



2e l 



d 3 p 
(2^)3 

d 3 p p\ 
(2tt) 3 E 3 



(f e (E)+f s (E)) 



1 

lu 2 E 2 - \e\ 2 B 2 



E 



2E 



(f e (E) + fc(E)) 



ir P {uj) = 2e 2 to(\e\B) 



d 3 p 
(2^p 



(f e (E) - f- e {E))- 



1 



2 E 2 -\e\ 2 B 2 



1-H 

E 2 



(5.33) 



where we can replace p 2 — > 2p 2 /3 and E — > \p\ in the integrand. The angular integration is trivial but the remaining 
integration over p cannot be carried out exactly in general. We therefore consider the classical and the degenerate 
limits, separately. 



1. Case I - Classical Gas 



In this case 

\ > M , (5-34) 

and 

f e , s (E) = e-ft**") . (5.35) 
The relationship between the chemical potential and the electron charge density is then given by 

n e - n s = —^—r sinh/3/x ~ , (5.36) 

7r z p J 7r z p^ 

and in addition 

4 4 

n e + n g = -—. rCOsh/3/i ~ — — . (5.37) 
The integration over p is also trivial for ttl and yields 

e 2 /3 

ttl = — (n e + ng) . (5.38) 

On the other hand, for ttt.p the integration over p cannot be done exactly. The formulas for these quantities can be 
expressed in the form 

ttt = -T- K + n s )F 3 

3 \ u> 

e 2 p n B (n B 

hp = — (n e -n g ) F 2 [ , (5.39) 



LO V LU 
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where 

F n (z) = / dxe~ x - , (5.40) 

Jo % - z 2 

with 

n B = \e\B0. (5.41) 

The quantity plays in the present case the role that ljb occupies in the non-relativistic case. We can obtain 
approximate values for the functions F n (z) by considering the region in which z > 1 or z C 1. 

The term e~ x x n is dominated by the values of x f» 1. For z 3> 1, the singularity at x — z then lies outside the 
effective region of integration, and we can substitute 

' ' (5-42) 



2 ' 



Therefore 



1 - z 



Fniz)-^—^, (5.43) 



and, from Eq. (5.39), we can write 



7r T = 2e (3{n e + n g ) 

,2 



uj 2 -n 2 



ttp = -£(n e - rig) -5— It > ( 5 - 44 ) 



3 v c c/ Vw 2 -^ 

for fis ^> w. In the opposite regime z« 1, the singularity is within the effective region of integration and requires 
some care for finite z. In the extreme limit, 

F„(0) = (n-2)! (5.45) 

and 

TTT = —\ n e + n e) > 

np = — {n e -n s ) , (5.46) 

0J 

for FIb <C ui. These expressions correspond to the weak- field limit used in Ref. jyj. 
For arbitrary values of z, the functions can be calculated by substituting the equalities 

x 2 i z ( 1 1 



2 \x — z x + z 
x 3 -,- / 1 I 



(5.47) 



x 2 — z 2 2 \x — z x + z / 

in Eq. ( p.4C| ), followed by the change of variable t — x — z. In this way, F 2 and F 3 are found to be given by 

F 2 (z) = l- 7 -[e- z Ei{z)-e*Ei{-z)\ , 
z 2 

F 3 {z) = 1-— [e- z Ei(z) + e z Ei(-z)\ , (5.48) 
in terms of the Exponential Integral function 

f* e* 

Ei{z) = V I —dt. (5.49) 

•J — oo t 
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By means of Eq. J53g| ), 7Tt,p can then be computed for any value of z. In particular, using the expansion 



^(z) = ln(|z|) + 7 + z+ T + - + 0(z 4 ) 

where 7 = 0.577215 is the Euler constant, for z small we obtain 

F 2 (z) = 1- z 2 + z 2 ( 7 + ln|z|) + 0(z 4 ), 
F 3 (z) = l-z 2 (ln(|z|)+ 7 ) + 0(z 4 ). 



(5.50) 



(5.51) 



When these expressions are substituted in Eq. fl5.39| ) they of course reproduce the lowest order terms given in Eq. 
( |5.46 ) plus the 0{Vt 2 B ) corrections to them. 



2. Case II - Degenerate gas 



This case corresponds to the limit 



The distribution functions are 



(5.52) 



f e = ©(pf -P) , 

h = 0, 



(5.53) 



where the Fermi momentum is the same as the chemical potential, and its relationship with the electron number 
density is 



The formula for ttl in Eq. ( p.33|) yields 



while the formulas for ~kt.p become 



7TL 



3tt 2 



e 2 n e 

PF 



(5.54) 



(5.55) 



e 2 \ 2 



PF 



dp- 



P 



7T 2 7 3 Jo p 

e\B\ 1 



2 _ e 2 B 2 ' 



UI J 3 



dp- 

p 



2 _ e 2 B 2 ' 



(5.56) 



With the instruction that the integrals in Eq. (5.56) are to be evaluated by taking the principal value part, they yield 

uj 2 

- log 1 - - 



e 2 n e 



Hp 



PF 

e 2 n e \ flp 



1 ^Fi 

1 + -r log 



i + ^iog 



n, 



where 



(5.57) 



(5.58) 



20 



Weak-field limit 



By this we mean that yfejl? is sufficien tly sm all, c ompa red to an average electron energy (E e ) and u>, so that we 
can literally take the B — > limit in Eqs. (4.30) and (4.31), and keep only the term that is linear in B. In this case, 
Eqs. Q5.271 ), Q5.28D and ( ^29| ) hold. As shown in Appendix g, using them in Eq. ( [4.30[ ) and taking the B — » limit 
results in 



ir T (uo) = 4|e| 

ir P {u) = 2|e| 2 

TTifw) = 4|e| 2 

2^2 



d 3 v E 2 - ±n 



(2tt) 2 2E 
2\e\B\ f d 3 p 



E 2 -u 2 /4 ' 
E 



3E 



(2K) 2 2E {h{E) - fs{E)) E 2 -u 2 /A 



d 3 n E 2 



(2n) 2 2E' 



E 2 -uj 2 /A ' 



where we have put p 2 — > |p in the integrand. For w <C 2(E e ), these reduce further to 



7t t (cj) = 7r L (o;) = = 4|e|" 



d 3 p 



-(f e (E)+f g (Ef)) 1 



(2tt) 2 2S 

which is the standard expression for the plasma frequency squared, and 

'4|e| 3 S\ /■ d 3 p {f e {E)-f- e {E)) 



P 

ZE 2 



7Tp(w) 



(2tt) 2 2£; 



2p 2 
3^ 2 



(5.59) 



(5.60) 



(5.61) 



both of which can be evaluated explicitly for the various limiting forms of the electron distribution functions that 
we have considered. We do not proceed any further in this direction, but we mention that the results thus obtained 
agree with the corresponding ones obtained in the previous sections, in th eir co mmon range of validity. That is, for 
example, for a degenerate distribution in the relativistic regime, from Eq. (5.61) we obtain 



\e\ 3 B PF 
3tt 2 uj 



(5.62) 



which agrees with the formula that is obtained from Eq. (5.57) in the limit flF <C w (which corresponds to the 
weak- field limit there ). For a non-relativistic gas, putting f s w and neglecting the term p 2 /E 2 in the integrand of 
Eqs. (ET601) and fcM) , 



irp(uj) 



ir T (u>) , 



(5.63) 



with wp = which coincides with the result obtained by taking the weak-field limit (uip <§C to) in Eq. (5.10). 
Analogous relations can be verified for other cases as well. However, we e mphasize that the results of the previous 
sections hold for a wider range of conditions than those for which Eq. (5.59), and the relations based on it, hold. The 
latter are limited to situations in which retaining only the linear terms in B is a valid approximation. 



D. Faraday Rotation 



The fact that the two transverse modes have different dispersion relations leads to the Faraday rotation effect, as 
is well known. After traveling a distance L, the direction of polarization of the wave has rotated by an angle 



= — LoAnL , 



(5.64) 



where 



An = (n_(w) — n + (w)) , 



(5.65) 



with n^(w) being the refractive indices of the left and right polarized modes, respectively. These can be computed 
from the dispersion relations given above by using n(uj) = Q/lo. For illustrative purposes, taking as an example the 
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case of a non-relativistic gas, for propagation parallel to the magnetic field from Eq. (5.14) we find 

0* 



n 2 » = 1- ° 

UJ{U) — U>b) 

n 2 _(uj) = 1- 



O 2 

"0 



and therefore 



Writing 



n_(ui) — n+(uj) 



ui(ui + uj b ) ' 



An = 



we then obtain 



n 2 _(uj) — n 2 + {ui) 
n_(w) + ' 

L 2n(w)(LU 2 -w 2 B )' 

where we have defined n(tu) — + n +(u>) )/2. 

For values of u> such that w 3> Oq, Eq. (5.66) tells us that n_ ~ n+ ~ 1 and therefore 



L 2(u 2 -oj 2 b ) 



(5.66) 



(5.67) 



(5.68) 



(5.69) 



(5.70) 



This result coincides with the one given in Ref.|16|. On the other hand, for w > wg, which corresponds to the 
weak- field approximation [p. 8| , 



9 £1qU^ 
L ~ 2uuj 2 



where, from Eq. (5.66), 



(5.71) 



(5.72) 



For other values of w, the formula in Eq. (5.69) interpolates nicely between the two limiting cases we have mentioned. 

More gener ally, for propagation parallel to the magnetic field as we have considered above, the rotation angle is 
given by Eq. (5.64) with 

_ ,_ 1 ±p L 3 

Thus, explicit formulas for the angle of rotation can be deduced simply for the other situations of interest using the 
results for 7Tt,p that we have already obtained for the various cases. 



VI. DISCUSSION AND CONCLUSIONS 



The subject of the propagation of a photon in a background of particles embedded in an external magnetic field 
appears in many physical contexts. As we have mentioned in the Introduction, there exist situations of interest for 
which neither the semiclassical methods nor the linear field-theoretic approach are directly applicable. In these cases, 
a more general field-theoretical treatment that does not involve the weak field assumption is required. 

With this motivation in mind, in this work we have reexamined the subject. We have given a general decomposition 
of the photon self-energy in a matter background that contains a magnetic field, in terms of the minimal set of tensors 
consistent with isotropy and the transversality condition. From this result, we have shown that the self-energy can be 
expressed in terms of nine independent form factors, that in the long wavelength limit reduce to three. In this limit, 
by applying the (real-time) finite temperature field theory method, we have calculated the one-loop formulas for the 
form factors. The formulas obtained in this way are valid for arbitrary distributions of the electron gas and for strong 
magnetic fields. They were explicitly evaluated for a variety of conditions of interest in physical applications, including 
the weak-field limit as a special case. From them we determined the photon dispersion relations and computed the 
Faraday rotation for various cases. They reproduce the well-known semi-classical results when the appropriate limits 
are taken, but they remain valid for more general situations, including those in which the magnetic field is not weak. 



22 



APPENDIX A: THERMAL ELECTRON PROPAGATOR 



In this section we explain briefly the formula for the electron propagator that will be used in the subsequent 
calculations. Our staring point is the expression given in the book by Itzykson and Zuberjlj]]. In vacuum, but in the 
presence of a constant magnetic field, the electron propagator in coordinate space is given by 



Sa(x, x') = [i$ x — eA'ix) + m](— i) / cItU(x,x';t) 



(Al) 



where [n9 



U (x, x'; t) 



16tt 2 t 2 
1 



')exp \ ~(x- x')eF coth(eFr)(x -x')+i Qer ■ F 



+ m — le I t 



Tr In [(eFr)" 1 sinh(eFr)] 



and with the convention that e is the charge of the electron. The factor <f>(x, x') is given in general by 



(x, x) — exp |— ie J d£ 



A^) + \f^-x' 



and it depends on the gauge but is independent of the path of integration, and it has the form 



(x,x')=exp\-F^x ,l/ \ , 



in the gauge in which 



A„ 



The expression for the propagator in momentum space is obtained by using the formula 



e 4 



TrVDet(A) ) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



where T)et(A) stands for the determinant of the matrix formed by the coefficients „. For the ensuing manipulations 
it is useful to observe that 



(A7) 



where P M „ is defined in Eq. ( 2.15 ). Using the multiplication rules P 2 = R and R 2 = R, the following relations are 
readily derived, 



cosh(eFr) = 1 + R [cos(s) - 1] 

1 -1 

- 1 



[cosh(eFr)] 1 = I + R 
{eFry 1 sinh(eFr) = 1 + R 



cos(s) 
sin(s) 



(A8) 



from which the rest of the formulas that we need can be obtained. In particular, 

tan(s) 



[(eFr) coth(eFr)] " 1 = 1 + R 
Det([(eFr)coth(eFr)]- lA - ^ tan(s) 



Trln[(eFT)" 1 sinh(eFr)] = In 



sin(s) 



(A9) 
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where s is the varia ble d efined i n Eg . (|3.6|) 



Thus, using Eqs. ( [Aq ) in Eq. ( |A2|) , with the help of Eq. dA9| ) we obtain 



U{x, x'; t) = 4> 



1 



d 4 p 



cos(s) J (2tt) 



e -ip-(x-x') e ~iT[<!>(p,s)-% a-F+ie] 



(A10) 



with $ as defined in Eq. (3.4) in the text. For the gauge choice given in Eq. (A5), which implies Eq. (A4), it follows 
that 



(i$ x - + m)U(x, x'\t) = <j) 



Using the relation 



which follows by noticing that 



together with the fact that (75 $5) 2 



where 



with 



1 



" —ip-(x-x) 



cos(s) J (2tt) 4 



\j> - tan(s)iP Ml ,7'y y + m] e 



-it \<b{p,s) — ftr-F+ic 



cos(e-Br) + 175^ sin(e£?r) , 



icr • F = 575^ 



1, Eq. ( |A1| ) then yields 
d 4 p 



SU(x,a;') = (j) 



(2tt)< 



POO 

tS F (p) = / drG{p, sy^'- 
Jo 



G(p, s) = \$- ta,n(s)iP fl ^p u + m] [1 + i7 5 #tan(s)] . 
In writing Eq. (A15), we have made the change of variable r — > — r. Using the identity 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 



G can be expressed in the final form given in Eq. ( |3.4| ) in the text. Finally, the 11 component of the thermal propagator 
is obtained from the formula 

S e (p) = S F (p) - (S F (p) - S F (p)) Ve , (A18) 

where r/ e is the background-depe nden t factor defined in Eq. ( |3.9| ) and, as usual, S F = j°S F ~f Q . This leads to the 
decomposition of S e given in Eq. (3.2), with the thermal part SVe as defined in Eq. (3.8). 

We would like to mention the following point. Strictly, the Fourier expansion given in Eq. (A10) holds only for those 
values of r for which cos(s) ^ 0. A t the points where cos(s) = 0, the expansion given in EqT^SlO) is not valid, and 
the inte grati on over r in Eq. ( A15 ) is not defined at those points as it stands. Therefore, we define the r integration 
in Eq. (A15), by deforming the line of integration to lie just below the real axis, so that the points in question are 
avoided. It is not difficult to see that this choice, as well any other that avoids those points (such as, for example, 
choosing the integration line to lie just above the real axis), amounts to a redefinition of the function 4>(x, x'). 



APPENDIX B: EVALUATION OF J 



In order to evaluate the integral Jo defined in Eq. (3.21), we first observe that, as a function of ip, Jq is analytic in 



the upper-half plane. Therefore, if we can evaluate it for a particular subregion of that plane, the result is valid for 
the entire region by analytic continuation. For the subregion Re ip < (or, equivalently, Euclidean momenta pu) we 
proceed as folows. 
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The integrand is an analytic function of s in the lower-half s-plane. Therefore we can deform the integration path 
and integrate over the negative imaginary axis. Then setting 

s = -it , (Bl) 

dte**e- a e^ , (B2) 

u = er 2t , (B3) 
2u 

itan(-ii) = 1. (B4) 



we have 



where 



e LB 



and we have used 



Using the generating formula of the Legerre polynomials 

xz oo 

n=0 



(valid for < x < oo and \z\ < 1), equation (B2) becomes 



Jo = 



e LB 



(B5) 



(B6) 



n=0 



where D n (a) is the function defined in Eq. ( [3.26| ) in the text. Equation (B6) is easily integrated term by term to yield 
the final result quoted in Eq. ( 3.25 ), which also holds for Re^ > (or equivalently, Minkowskian values of p||) by 
analytic continuation. 



APPENDIX C: LOW FREQUENCY LIMIT 



To arrive at Eqs. (5.4) and ( |5.5| ) we apply the low frequency conditions to Eq. (4.31) as follows. Let us define 

1 1 



± 



By combining the denominators, 



uj 2 + 2uE n - 2\e\B uj 2 - 2ujE n - 2\e\B 
2{uj 2 - 2|e|B) 



D+ = 
D_ = 



Now define the combinations 



A 



(±) 



{uj 2 - 2\e\B) 2 - Auj 2 E 2 ' 

-4uE n 

(u> 2 - 2\e\B) 2 - Auj 2 E 2 ' 



D+ ± (B -> -B) , 
Z?_ ± (B -> -B) , 



in terms of which the formulas for 7r^p given in Eq. (4.31) are expressed as 



uj) 



dp\\ 
E n 



4 0) m 



(f e (E n ) + f g (E n )) [ul 
dp„(f e (E n ) + f e (E n )) {ujD_} , 



L +) +2n|e|BA ( F +) } 



-T„ "'{Uj) 



|l(/ e (B„)-/ g (B„)){ 



uKA*;' + 2n\e\BA { ~ } 



*|e|BA^} , 



4 0) M 



dp\\(f e (E n )-f e (E n )){uD+} 



(CI) 



(C2) 



(C3) 



(C4) 
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From their definitions, and using Eq. (CI), 



2(w 2 - 2\e\B) 



± 



2(uj 2 + 2\e\B) 



(w 2 - 2\e\B) 2 - Auj 2 E 2 (uj 2 + 2\e\B) 2 - Auj 2 E 2 



A 



(±) 



-4u)E„ 



± 



(w 2 - 2\e\B) 2 - Alo 2 E 2 {uj 2 + 2\e\B) 2 - Auj 2 E 2 



(C5) 



It is seen that all the A's can be approximated by letting (to 2 — 2\e\B) 2 — > 4\e\ 2 B 2 in the denominators, exce pt A*i_ \ 
which would give zero. This in particular means that we cannot take the low frequency limit in Eq. ( 4.31 ) by just 
setting lu 2 — > in the denominators, or otherwise we would miss some important terms. In order to make a systematic 
expansion of A_ for uj <C ^/|e|B, we use the relation 



which follows from using the trivial identity 

1 1 



1 L0 2 -2\e\B 



loE„ 



2ujE„ 



D 



+ ■ 



x + y x-y y 



2-x 



Then going back to the definition in Eq. 



, we can re-express A_ as 



x+y x — y 

(-) 



A (-) = _( u 



\2E n 



\e\B_ 

UjE n 



,(+) 



(C6) 



(C7) 



(C8) 



Up to this point the relations are exact. Now we can proceed to make the approximation uj -C y|e]S systematically, 
by letting (uj 2 ± 2\e\B) -> ±2\e\B in Eq. (p|, except for A ( f ) , for which we use Eq. (p|. Thus, 



Similarly, from Eq. (C2) 



Ai +) 



At' = 



A<_ +) = 



A« = 



uj 2 E 2 - \e\ 2 B 2 

2\e\B 
uj 2 E 2 - \e\ 2 B 2 

2wE n 
cu 2 E 2 - \e\ 2 B 2 

2uj\e\B" 



1 



uj 2 E 2 - \e\ 2 B 2 



(C9) 



Da 



\e\B) 



uj 2 E 2 - \e\ 2 B 2 

UjE n 



2 E 2 



l B 2 



(CIO) 



The formulas quoted in Eq. ( p.4[ ) follow by substituting these approximate expressions for the A's and the D's into 
Eq. @. 



For tt^ 1 \uj) we assume that ui <C 2(E e ) to approximate 

1 



1 



uj 2 + 2ujE n uj 2 - 2ujE„ 



2EI 



(Cll) 



from which the formula in Eq. (5.5) follows. 
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APPENDIX D: THE WEAK-FIELD (LINEAR) LIMIT 



In order to obtain the weak-field limit formulas quoted in Eq. (5.59), we use the expressions given in Eq. ( C4j ) for 
n T ,p{u). From Eq. @, as B -> 



>(+) 



2ujE 



lu 2 (E 2 -lo 2 /A) ' 



Then, remembering that Eq. ( 5.29 ) holds, 



2|e| 



A 



(+) 



E 2 - lu 2 /A 



(2tt) 2 2E 



(f e (E)+U(E)) 



E 2 -lo 2 /A 



[2E 2 -pl] , 



(Dl) 



(D2) 



which yields th e fo rmula given in Eq. ( [5.59| ), after replacing p\ 
tt p . From Eq. 



A 



(-) 



Ids 



LJ 4 - Alo 2 E 2 

,2 



V2w J By + 



|p 2 in the integrand. We proceed similarly for 



Ids 



Ids 



8w 2 



w£ - AE 2 l 



where we have used Eq. (|CS|) for AL } . Then substituting in Eq. dC4|) 

'2|dS 



7rp(w) 



2d 



dP (fe(E) - fe(E))- E 



{2n) 2 2E yJCK ' JCK "E 2 -lo 2 /A 
Finally, for itl{u) we apply Eq. (5.29) directly in Eqs. ( 4.3CQ and (4.31) and obtain 



7TL (W) 



4d 



d 3 p 



{2it) 2 2E 



(f e (E)+f s (E))[4m 2 e + 4pj 



1 



1 AE 2 - to 2 



(D3) 



(D4) 



(D5) 



which can be expressed as quoted in Eq. (5.59). 



M. Giovannini, Phys. Rev. D 56, 3198 (1997). 

M. M. Giovannini and M. Shaposhnikov, Phys. Rev. D 57, 2186 (1998). 
M. J. Rees, Quart. J. R. A. S. 28, 197 (1987). 
P. P. Kronberg, Rep. Prog. Phys. 57, 325 (1994). 
T. Kolatt, Phys. Rev. D 56, 3198 (1997). 
A. Loeb and A. Kosowsky, Astrophys. J. 469, 1 (1996). 

G. Fleishman, Q. Fu, M. Wang, G.-L. Huang, and V. Melnikov, Phys. Rev. Lett. 88, 251101 (2002). 
T. Verbiest, M. Karunanen, and A. Persoons, Phys. Rev. Lett. 82, 3601 (1999). 
S. Ishimam, Statistical Plasma Physics (Addison Wesley, Massachusetts, 1992). 
E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics (Pergamon Press, Oxford, 1981). 
A. K. Ganguly, S. Konar, and P. B. Pal, Phys. Rev. D 60, 105014 (1999). 
K. H. et. al., Astrophys. J. Lett. 510, Llll (1998). 

E. V. Gotthelf, G. Vasisht, and T. Dotani, Astrophys. J. Lett. 522, L49 (1999). 
R. Narayan, P. Paczyhski, and T. Piran, Astrophys. J. Lett. 395, L83 (1992). 
T. Vachaspati, Phys. Lett. B 265, 258 (1991). 

M. Harwit, Astrophysical Concepts (Springer- Verlag, New York, 1998). 
C. Itzykson and J. B. Zuber, Quantum Field Theory (McGraw-Hill, New York, 1997). 

This result agrees with what can be deduced from the results given in Ref. [|ll| for the non-relativistic gas. Notice that the 
quantity denoted by d<b/dl in that reference corresponds to a; An in our notation, and is twice the rotation angle. 
[19] Our formula for U differs from the one quoted in the book cited, by a factor of e in the term ^cr • Ft, which we believe that 
it is missing there, and in the overall sign, which believe is an error in the determination of the normalization constant C 
in the book. 



